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Gravitational wave astronomy 

Gravity is one of the fundamental forces of Nature, and 
it is the dominant force in most astronomical systems. 
In common with all other phenomena, gravity must 
obey the principles of Special Relativity. In par- 
ticular, gravitational forces must not be transmitted or 
communicated faster than light. This means that when 
the gravitational field of an object changes, the changes 
ripple outwards through space and take a finite time to 
reach other objects. These ripples are called gravita- 
tional radiation or gravitational wavesl^ 

In Einstein's theory of gravitation (see General 
Relativity and Gravitation), as in many other 
modern theories of gravity (see Non-general Rel- 
ativity Theories of Gravity), gravitational waves 
travel at exactly the speed of light. Different theories 
make different predictions, however, about details, such 
as their strength and polarization. There is strong in- 
direct observational evidence (see Binary Stars as 
a Probe of General Relativity, Hulse- Taylor 
Pulsar) that gravitational waves follow the predictions 
of general relativity, and instruments now under con- 
struction are expected to make the first direct detec- 
tions of them in the first years of the 21st century. 

These instruments and plans for future instruments 
in space are described in the article Gravitational 
Radiation Detection on Earth and in Space. 
Detectors must look for gravitational radiation from 
astronomical systems, because it is not possible to gen- 
erate detectable levels of radiation in the laboratory. 
It follows that gravitational wave detection is also a 
branch of observational astronomy. 

The most striking aspect of gravitational waves is 
their weakness. A comparison with the energy in light 

*To be published in the Encyclopedia of Astronomy and Astro- 
physics (Institute of Physics Publishing, Bristol, and Macmillan 
Publishers Ltd, London, 2000). 

J They are also sometimes referred to as gravity waves, but 
since this term has a different meaning in meteorology and stel- 
lar hydrodynamics, we will avoid it here. See Solar Interior: 
Helioseismology. 



will illustrate this. The human eye has no trouble sens- 
ing the light from the planet Jupiter: the amount of en- 
ergy that passes through the iris of the eye is far more 
than the minimum the eye can detect. Yet several times 
a week a gravitational wave, generated in a far distant 
galaxy, carries a similar amount of energy into the eye, 
and we don't notice it. 

The reason is that gravity is the weakest of the fun- 
damental forces, and the disturbance created by even 
such an energetic wave is so tiny that no man-made in- 
strument has so far registered it. While all the energy 
in the light from Jupiter that enters the eye is absorbed 
in the eye, the gravitational wave passes right through, 
leaving behind almost none of its energy. All the mat- 
ter in the present Universe is similarly transparent to 
gravitational waves. 

Gravitational radiation is today one of the last un- 
opened windows into the Universe. There are at least 
five reasons motivating scientists to develop gravita- 
tional wave astronomy: 

• The weakness with which gravitational waves in- 
teract with matter is a great advantage for astron- 
omy. It means that gravitational waves arrive un- 
affected by any intervening matter they may have 
encountered since being generated. There is no 
significant scattering or absorption, although they 
will be deflected by a Gravitational Lens in 
the same way as light. Gravitational waves carry 
uncorrupted information even if they come from 
the most distant parts of the Universe or from its 
most hidden regions, like the interiors of Super- 

NOVAE. 

• Gravitational waves are emitted by the bulk mo- 
tions of their sources, not by individual atoms or 
electrons, as is normally the case for electromag- 
netic waves. They therefore carry a completely 
different kind of information about their sources 
from that which is normally available in observa- 
tions of Binary Stars, supernovae, and Neu- 
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tron Stars. 

• Gravitational waves can be emitted by black holes, 
which are described in the article General Rel- 
ativity and Gravitation. Indeed, gravita- 
tional waves provide only way to make direct ob- 
servations of these objects. Since there is now 
strong indirect evidence that giant black holes in- 
habit the centers of many (or even most) galax- 
ies (see Supermassive black holes in AGN), 
and since smaller ones are common in the Galaxy 
(see Black Hole Candidates in X-Ray Bi- 
naries), there is great interest in making direct 
observations of them. 

• Gravitational waves can come from extraordinar- 
ily early in the history of the Universe. The elec- 
tromagnetic radiation from the Big Bang is called 
the Cosmic Microwave Background. Obser- 
vations of it describe the Universe at it was about 
10 5 years after the Big Bang. Studies of cosmo- 
logical Nucleosynthesis give information about 
what the Universe was like as little as 3 minutes 
after the Big Bang. Gravitational waves, if they 
can be detected, would picture the Universe when 
it was only perhaps 1CP 24 seconds old, just at the 
end of Inflation. 

• Gravitational radiation is the last fundamental 
prediction of Einstein's general relativity that has 
not yet been directly verified. If another theory of 
gravity is correct, then differences could in prin- 
ciple show up in the properties of gravitational 
waves, such as their polarization. In principle, 
there must be a better theory of gravity, since 
general relativity is not a quantum theory, a defi- 
ciency that theoretical physicists today are work- 
ing hard to remedy. The majority belief today is 
that there should be a unified theory of the fun- 
damental forces, in which gravitation is related to 
the other forces. Evidence for the nature of this 
relation could show up in observations of grav- 
itational waves, particularly those from the Big 
Bang. 

These motivations and their implications are devel- 
oped in the following sections. Each section begins with 
an introduction to the physical ideas and then develops 
some of the mathematical details. 

The physics of gravitational radiation: 
weakness and strength 

The starting point for understanding gravitational ra- 
diation is Newtonian gravity. The weakness of gravity 
is evident. If a child picks up a book, she defeats the 
cumulative gravitational pull of the entire planet Earth 



on the book. The strength to do this comes from the 
chemical forces in her muscles, which come from elec- 
tromagnetic interactions. 

In fact, the electromagnetic force between the elec- 
tron and the proton in a hydrogen atom is 2 x 10 39 times 
bigger than the gravitational force between them. The 
reason that gravity can nevertheless dominate on the 
cosmic scale is that opposite electrical charges cancel 
each other, while the gravitational forces of all the par- 
ticles add. 

Another fact about gravity that was known to New- 
ton is what is now called the equivalence principle (see 
General Relativity and Gravitation). This is 
the principle that all bodies accelerate in the same way 
in a gravitational field, so that the trajectory that a 
freely falling body (a body influenced only by gravity) 
follows in a given gravitational field depends only on its 
starting position and velocity, not on what it is made 
of. 

Imagine now a machine made in some way to detect 
gravitational waves. Whatever the method of detec- 
tion, a wave needs somehow to alter the internal state 
of the detector. If the wave carries a gravitational field 
that is completely uniform across the detector, then by 
the equivalence principle all of the parts of the detector 
will accelerate together, and its state will not change at 
all. To detect a gravitational wave, the machine must 
measure the non-uniformities of the gravitational field 
across a detector. 

These non-uniformities are called tidal forces, be- 
cause they produce the stretching effects that raise tides 
on the Earth. Gravitational waves are traveling tidal 
forces. 

Newton's theory of gravity had no gravitational 
waves. For Newton, if a gravitational field changed in 
some way, that change took place instantaneously ev- 
erywhere in space. This is not a wave. Let us consider 
what we mean by the term "wave" in ordinary language. 

Imagine a child's rubber duck floating in a bath tub 
half full of water. If a child presses down on the duck 
very gently, until is is nearly submerged, then the level 
of the water will rise everywhere in a nearly uniform 
way, and this is not called a wave. If instead he drops 
the duck, then the disturbance rises around the base 
of the duck rapidly, moves away from it, and eventually 
reaches the walls. This is a wave. Wave motion requires 
a finite speed for the propagation of disturbances. If the 
disturbance is very slow, as for the floating duck, then 
the wavelength is very long, and near the site of the 
disturbance the wave motion is not noticeable. We say 
we are in the "near zone" . But when we are more than 
a wavelength away, then we see waves, and this is the 
"wave zone" or "far zone" . For the dropped duck, we see 
that waves because their wavelengths are shorter than 
the size of the bath. In general relativity, the speed 
of gravity is the speed of light. Because of this finite 
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speed, gravity must exhibit wave effects. 

Many of Newton's contemporaries were unhappy 
that his theory of gravity was based on instantaneous 
"action at a distance" , but Newton's theory fit the ob- 
servational facts. If gravity had a finite speed of propa- 
gation, there was no evidence for it in the solar system. 
Interestingly, the brilliant 18th century French math- 
ematician and physicist Laplace tried out a variation 
on Newton's theory in which gravity was represented 
by something "flowing" out of its source with a finite 
speed. He reasoned that a planet like the Earth, moving 
through this fluid of gravity, would experience friction 
and gradually spiral in towards the Sun. 

Laplace could show that the observational limits on 
this inspiral even in his day were so stringent that the 
speed of gravity in his model needed to be huge com- 
pared to the speed of light. He did not find this result 
attractive and took the theory no further. (Laplace also 
explored the notion of what we now call a black hole, 
which for him was a region where gravity was strong 
enough to trap light.) 

It is interesting that today, observations of the two 
neutron stars in the binary system PSR1913+16 spi- 
raling together as they orbit one another provides the 
most convincing evidence that gravitational waves ex- 
ist and are as described in general relativity. (See below 
and the article Binary Stars as a Probe of Gen- 
eral Relativity.) Laplace had the right effect, but 
the wrong theory. This evidence is described in the next 
section. 

In general relativity, Einstein used the principle of 
equivalence as the basis for a geometrical description of 
gravity. In the four-dimensional world of space-time, 
the trajectory of a particle falling freely in a gravita- 
tional field is a certain fixed curve. Its direction at any 
point depends on the velocity of the particle. The equiv- 
alence principle implies that there is a preferred set of 
curves in space-time: at any point, pick any direction, 
and there is a unique curve in that direction that will be 
the trajectory of any particle starting with that veloc- 
ity. These trajectories are thus properties of space-time 
itself. 

Moreover, if there were no gravitational field, the 
trajectories would be simple straight lines. Even in a 
gravitational field, a small freely falling particle does not 
"feel" any acceleration: its internal state is the same as 
if there were no gravity. Therefore Einstein postulated 
that a gravitational field made space-time curved, and 
that the preferred trajectories were locally straight lines 
that simply changed direction as they moved through 
the curved space-time, in much the same way as a great 
circle on a sphere changes direction relative to other 
great circles as one goes along it. For weak gravitational 
fields of slowly moving bodies, Einstein's theory reduces 
to Newton's in the first approximation. 

For gravitational waves, one could make a very sim- 



ple detector just by monitoring the distance between 
two nearby freely falling particles. If they are genuinely 
free, then any changes in their separations would indi- 
cate the passage of a gravitational wave. Because this 
measures a tidal effect, the bigger the separation of the 
particles, the bigger will be the change in their separa- 
tion, at least for particles that are separated by less than 
a gravitational wavelength. Most modern gravitational 
wave detectors are designed to be as big as cost and 
practicality allows. These are described in the article 
Gravitational Radiation Detection on Earth 
and in Space. 

Although gravitational radiation is well understood 
in theoretical terms in general relativity, the complex- 
ity and non-linearity of Einstein's equations means that 
calculations are often difficult. In the historical de- 
velopment of general relativity, between 1915 and the 
1950's and 1960's, these mathematical difficulties cre- 
ated confusion over the physical nature of gravitational 
radiation, and in particular over whether they carried 
energy away from the source. Improved mathematical 
techniques finally resolved the matter in favor of the 
simple physical picture presented here, but this picture 
would not be complete without the strong mathemati- 
cal underpinning that now exists. 

The question of energy in gravitational waves is still 
a delicate one. There is no question that waves carry 
energy (and momentum) away from their sources. Nev- 
ertheless, it is not possible in general relativity to local- 
ize the energy in the radiation to regions smaller than 
about a wavelength. Indeed the equivalence principle 
shows that "point" particles feel nothing, no matter 
how strong the wave. The wave only acts by stretching 
space-time, producing a tidal distortion in the separa- 
tions between particles (see the discussion of polariza- 
tion below). 

For this reason, energy is localized only in regions, 
not at points. It is nevertheless real energy: the nonlin- 
earity of general relativity allows waves to create grav- 
itation themselves. Recent numerical simulations have 
shown that focussed gravitational waves can actually 
form black holes, trapping themselves. If the waves are 
weak, they enter the focussing region and re-emerge. If 
they are strong enough, they enter and never leave. 

Gravitational waves in a quasi-Newtonian model 

It is possible to calculate the approximate size of the 
effect of a given gravitational wave by beginning with 
Newtonian gravity and adding waves to it. In Newto- 
nian gravity the gravitational field produced by a mass 
M at a distance r is given by 

<j> = -GM/r, (1) 

where G is Newton's gravitational constant. The field 
of a gravitational wave must be a ripple on this, which 
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means a small change that oscillates in space and time. 
A suitable form for a change that propagates at the 
speed of light in the z-direction with an angular fre- 
quency u> is: 



GM , 

-e sm|w(z/c — t)J, 



(2) 



where e is a dimensionless number that would be ex- 
pected to be small compared to 1. Its size is the subject 
of the next main section. 

The field 5<f> produces an acceleration in the z- 
direction that depends on its z-derivative. Both 1/r 
and the sin() term depend on z. The derivative of 1/r 
will be proportional to 1/r 2 , which is how the acceler- 
ation falls off in Newton's theory (where is the only 
field). But the derivative of the sin() term does not 
change the 1/r; rather, it essentially just multiplies 8<j> 
by lo/c. At sufficiently large distances from the source, 
this term will dominate the 1 /r 2 term and the acceler- 
ation produced by the wave will be: 



Sa z — ~eu cos[u>(z/c — t)]. 



(3) 



Note that this term would not be present in the x- and 
y-derivatives, so these components of the acceleration 
are much smaller in this quasi-Newtonian model of a 
gravitational wave. 

Effect on a simple detector 

The tidal part of this acceleration, for a detector that 
has size £ in the z-direction, is to a first approximation 



d 9 GlStl 

£ ■ ~rSa z = euj £ — ~- sin[u;(z/c - f)J 



dz 



rc* 



(4) 



If a detector consists of two freely falling particles with 
this relative acceleration, the equation of motion for 
their separation £ will be 

d2£ 2„GM .... 

^2 = eu) bih[w(z/c - t)], (5) 

The dimensionless coefficient e(GM/rc 2 ) is typically 
very small. Even if e is of order 1, the other number is, 
with reasonable values for M and R, 



GM n A 21 / M 

— — = 2.4 x 10 I — — 

IMra 



re 1 



20 Mpc 



The mass and distance scales here are those appropriate 
to a neutron star in the nearest large cluster of galax- 
ies, the Virgo Cluster. (The distance unit is based on 
the astronomers' parsec, denoted pc, which is about 
3 x 10 16 m. The unit Mpc is a megaparsec.) It is be- 
lieved that several neutron stars are formed in the Virgo 
Cluster each year in supernova explosions. Ever since 
the beginning of the development of gravitational wave 



detectors such events have been high on the list of pos- 
sible sources of gravitational waves. 

To solve Equation (||), one takes r and z as con- 
stants, and uses the smallness of the right-hand-side, 
which implies that the changes in £ are tiny compared 
to t itself. On the right-hand side one can therefore 
replace £ by £q, the initial value of £, and then sim- 
ply integrate twice in time to get (for an initial value 
di/dt = 0) 



e(t) - £ GM 

7 = -e — T sm N z / c ~ *) ■ 

£o r& 



(6) 



The right-hand-side of Equation (Q) is identical to that 
of Equation (||) . This is an important conclusion which 
fits neatly with Einstein's geometrical conception of 
gravity: the size of a gravitational wave gives directly 
the stretching of the distance between nearby free parti- 
cles. It is conventional to call this h/2 and refer to h as 
the gravitational wave potential 



*:=2^ = 2f*. 
£o 

The amplitude of the oscillations of h is 



2c 



GM 



r<r 



(7) 



(8) 



It is evident from this that a detector must be able 
to measure changes in its own size that are smaller than 
one part in 10 21 to have a reasonable chance of making 
astronomical observations. The extraordinary small- 
ness of this effect also explains why ordinary objects 
in the Universe are transparent to gravitational waves. 
As the waves pass through them, they disturb them so 
little (parts per 10 21 typically) that the transfer of en- 
ergy to the object and any back-reaction effects of this 
on the wave are negligibly small. 

Energy flux carried by waves 

The energy in the waves can also be estimated from 
these equations and general physical principles. Quite 
generally, in classical field theories, the energy flux of 
a propagating sinusoidal plane wave is proportional to 
the square of the time-derivative of the fundamental 
field. In electromagnetism, the Poynting flux is propor- 
tional to the square of the time-derivative of the vector 
potential. 

In general relativity, the flux is therefore propor- 
tional to the square of the time-derivative of hit). The 
proportionality constant must be built only out of c, G, 
and pure numbers. To get the right units, it must be 
proportional to c 3 /G; to get the pure number, a calcula- 
tion in general relativity is required: 1 /32n for a linearly 
polarized wave. (Polarization is described later.) This 
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(9) 



/ 



100 Hz 



10 



-22 



W m~ 



for a wave with frequency / = lu/2tt. For comparison, 
reflected sunlight from Jupiter has a flux on Earth of 
2.3 x 10" 7 W m~ 2 , almost 100 times smaller than that 
of a gravitational wave with an amplitude of 10~ 22 ! 

Deficiencies of the quasi-Newtonian model 

The calculation and equations in this section have been 
framed within a modified Newtonian model of gravity 
with a propagation speed of c, and one would expect 
some differences from general relativity. The most im- 
portant difference is in the direction in which the tidal 
forces act. In the simple model, wave accelerations act 
in the z-direction, which was the direction of propaga- 
tion of the wave. This is called a longitudinal wave. 

In general relativity, gravitational waves are trans- 
verse waves: if the wave propagates in the z-direction 
then the tidal forces act only in the iry-plane. We will 
discuss later the exact form that their action in this 
plane takes. Remarkably, the rest of the formulas above 
are good approximations even in general relativity, pro- 
vided e is calculated correctly, as described in the next 
section. 

The emission of gravitational waves 

The previous section described the propagation of gravi- 
tational waves, their interaction with detectors, and the 
energy they carry. This section deals with the strength 
with which waves are emitted by astronomical bodies. 

In Newtonian gravity there is a fundamental the- 
orem, proved by Newton, that the gravitational field 
outside a spherical body is not only spherical, but it is 
the same as that of a point mass located at the origin of 
the body. It has the form given in Equation (|Tj). In par- 
ticular, the field is independent of the size of the body, 
as long as we consider only points outside it. This is 
true even if the star pulsates in a spherical manner. 

This theorem is essentially the same in general rel- 
ativity, and is known as Birkhoff's Theorem. Outside a 
spherical body the field is the same as that of a black 
hole of the same mass as the body (the Schwarzschild 
metric), even if the body is pulsating spherically. But 
if the pulsation is nonspherical, then the outside field 
will change. In general relativity the changes generally 
propagate as a wave. So gravitational waves will be 
emitted by nonspherical motions. 

In general the calculation of the emitted waves is 
extremely difficult, since the field equations of general 



relativity are a system of many coupled, nonlinear, par- 
tial differential equations. But in four circumstances the 
emission mechanisms are understood in some detail: 

• Small-amplitude pulsations of relativistic stars 
and black holes. Normally gravitational radiation 
carries away energy and damps pulsations away, 
but in rotating stars the opposite may happen: 
the radiated loss of angular momentum may allow 
the star to spin down to an energetically more fa- 
vored state, in which case the perturbation will 
grow, at least until nonlinear effects intervene. 
Discovered by S. Chandrasekhar and now called 
the Chandrasekhar- Friedman- Schutz (CFS) insta- 
bility, it is thought to limit the rotation speed of 
young neutron stars (see below) . Black holes also 
emit gravitational radiation when they are dis- 
turbed, e.g. by something falling into them, but 
they are not unstable: they always settle down 
into a steady state again. 

• Radiation from "test" objects orbiting black 
holes. If the mass of the object is small enough 
then the total gravitational field may be treated 
as a linear perturbation of the exactly known field 
of a black hole (the Schwarzschild or, with rota- 
tion, the Kerr solution). These studies give in- 
sight into the general problem of gravitational ra- 
diation, and they also predict gravitational wave- 
forms that might be observed by space-based ob- 
servatories looking at compact stars falling into 
the giant black holes in the centers of galaxies. 
(See below and the article Gravitational Ra- 
diation Detection on Earth and in Space.) 

• Weak gravitational fields and slow motion. Such 
weakly relativistic sources are studied in the post- 
Newtonian approximation, which includes higher- 
order corrections to Newtonian gravity from gen- 
eral relativity. This is analogous to the slow- 
motion multipole approximation that is so pow- 
erful in the study of electromagnetic radiation. 
Most realistic gravitational-wave sources can be 
studied to some approximation this way. 

• Collisions of black holes and neutron stars. These 
events, which are expected to be observed by 
gravitational wave detectors (see below) , must be 
modeled by solving the full set of Einstein equa- 
tions on a powerful computer. Techniques to do 
this are advancing rapidly, and simulations of re- 
alistic mergers of stars and black holes from in- 
spiraling orbits can be expected to yield useful 
results in the first years of the 21st century. 
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Quadrupole approximation 

The post-Newtonian approximation has so far been the 
most powerful of these methods, and it yields the most 
insight into the emission mechanisms. Its fundamental 
result is the quadrupole formula, which gives the first 
approximation to the radiation emitted by a weakly rel- 
ativistic system. 

The quadrupole formula is analogous to the dipolc 
formula of electromagnetism. In this language, 
monopole means spherical, which emits no radiation. 
This is also true in electromagnetism, where it is linked 
to conservation of charge. The "monopole moment" in 
electromagnetism is the total charge of a system, and 
since that does not change, there can be no spherical 
radiation. 

Again in electromagnetism, the dipole moment is 
defined as the integral 



pXid x, 



where p is the charge density and Cartesian co- 

ordinate. If this integral is time-dependent, then the 
amplitude of the electromagnetic waves will be propor- 
tional to its first time-derivative ddjjdt, and the radiated 
energy will be proportional (as we remarked earlier) to 
the square of the time derivative of this amplitude, i. e. 

to£J^M 2 l 2 - 

In the post-Newtonian approximation to general rel- 
ativity, the calculation goes remarkably similarly. The 
monopole moment is now the total mass-energy, which 
is the dominant source of the gravitational field for 
non-relativistic bodies, and which is constant as long 
as the radiation is weak. (Radiation will carry away 
energy, but in the post-Newtonian approximation that 
is a higher-order effect.) The dipole moment is given 
by the same equation as above, but with p interpreted 
as the density of mass-energy. 

However, here general relativity departs from elec- 
tromagnetism. The time-derivative of the dipole mo- 
ment is, since the mass-energy is conserved, just the 
integral of the velocity Vi : 



pvid x. 



(10) 



But this is the total momentum in the system, and (to 
lowest order) this is constant. Therefore, there is no en- 
ergy radiated due to dipole effects in general relativity. 
The gravitational field far from the source does con- 
tain a dipole piece if di is non-zero, but this is constant 
because it reflects the fact that the source has non-zero 
total momentum and is therefore moving through space. 

To find genuine radiation in general relativity one 
must go one step beyond the dipole approximation to 
the quadrupole terms. These are also studied in elec- 
tromagnetism, and the analogy with relativity again 



is close. The fundamental quantity is the spatial ten- 
sor (matrix) Qjk, the second moment of the mass (or 
charge) distribution: 



Qjk = J pxjX k d 3 x. 



(11) 



A gravitational wave in general relativity is represented 
by a matrix hjk rather than a single scalar h, and its 
source (in the quadrupole approximation) is Qjk- 

As in electromagnetism, the amplitude of the radi- 
ation is proportional to the second time-derivative of 
Qjk, and it falls off inversely with the distance r from 
the source. A factor of G/c 4 is needed in order to get a 
dimensionless amplitude h, and a factor of 2 to be con- 
sistent with the definition in Equation (||). The result 
for hjk is: 

2Gd 2 Q 3k 



hjk - 



(12) 



rc 4 dt 2 

General relativity describes waves with a matrix be- 
cause gravity is geometry, and the effects of gravity are 
represented by the stretching of space-time. This ma- 
trix contains that distortion information. Here is the 
information about the transverse action of the waves 
that the quasi-Newtonian model of the last section did 
not get right. 

Simple estimates 

If the motion inside the source is highly non-spherical, 
then a typical component of d 2 Qj k /dt 2 will (from Equa- 



tion (|11|)) have magnitude Mv^ s , where v 2 N s is the 
non-spherical part of the squared velocity inside the 
source. So one way of approximating any component 
of Equation (112) is 



2GMv 2 NS 
re 4 



(13) 



Comparing this with Equation (||) we see that the ratio 
e of the wave to the Newtonian potential is simply 



J N.S. 



By the virial theorem for self-gravitating bodies, this 
will not be larger than 



e < 4>int/c 2 , 



(14) 



where <^ n t is the maximum value of the Newtonian 
gravitational potential inside the system. This provides 
a convenient bound in practice. It should not be taken 
to be more accurate than that. 

For a neutron star source one has <pi n t ~ 0.2c 2 . If 
the star is in the Virgo cluster, then the upper limit on 
the amplitude of the radiation from such a source is 5 x 
10 -22 . This has been the goal of detector development 
for decades, to make detectors that can observe waves 
at or below an amplitude of 10 -21 . 
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Polarization of gravitational waves 

The matrix nature of the wave amplitude comes from 
general relativity and has no Newtonian analog. In or- 
der to find the effect of the waves on the separation 
of two free particles (the idealized detector), one has 
to start with hjk as given by Equation (|l2j ) or by any 
other calculation, and then do three things: 

1. Project the matrix hjk onto a plane perpendic- 
ular to the direction of travel of the wave. In 
the simple case considered above, where the wave 
was traveling in the z-direction, this means leav- 
ing the components {h xx , h xy , h yy } alone and 
setting the remaining components to zero. It is 
then a two-dimensional matrix in the transverse 
plane. 

2. Remove the two-dimensional trace of the pro- 
jected matrix. Call the resulting matrix hJjF , 
where TT stands for Transverse- Traceless. In 
the example this means subtracting (h xx + h yy )/2 
from both h xx and h yy . Then there are only two 



independent components left, h 



'TT 



-h TT 



TT 
xy 



■TT 



and 



3. To find the change in the separation of two par- 
ticles that have an initial separation given by the 



vector let the matrix hjk act on it 



SI 



(15) 



It is clear that any longitudinal component of the 
separation lj between the particles is unaffected by the 
wave (in the example, this is the z-separation) , and that 
there are two degrees of freedom (the two independent 
components of hJjF) to move particles in the plane per- 
pendicular to the propagation direction. These two de- 
grees of freedom are the two polarizations of the wave. 

Fig. |l| shows the conventional definition of the two 
independent polarizations, from which any other can be 
made by superposition. What is shown is the effect of 
a wave on a ring of free particles in a plane transverse 
to the wave. The first line shows a wave with h xy = 0, 
conventionally called the "+" polarization. The bottom 
line shows a wave with h xx — 0, the "x" polarization. 

Luminosity in gravitational waves 

The energy carried by the gravitational wave must be 
proportional to the square of the time-derivative of the 
wave amplitude, so it will depend on the sum of the 
squares of the components d 3 Qjk/dt 3 . The energy flux 
falls off as 1/r 2 , but when integrated over a sphere of 
radius r to get the total luminosity, the dependence on 
r goes away, as it should. The luminosity contains a fac- 
tor G/c 5 on dimensional grounds, and a further factor 
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Figure 1: Polarization of gravitational waves. The cen- 
ter line gives the wave as a function of time, with an 
amplitude of h = 0.2, and the top and bottom lines 
show to scale the distortions produced by two polariza- 
tions with this amplitude. 

of 1 /5 comes from a careful calculation in general rela- 
tivity. The result is the gravitational wave luminosity 
in the quadrupole approximation: 



L ~ — 

gw ~ 5c 5 



(16) 



where Q is the trace of the matrix Qjk- Its squared 
third derivative must be subtracted in order to ensure 
that spherical motions do not radiate. 

This equation will be used in the next section to 
estimate the back-reaction effect on a system that emits 
gravitational radiation. 

Emission estimates 

Until observations of gravitational waves are success- 
fully made, one can only make intelligent guesses about 
most of the sources that will be seen. There are many 
that could be strong enough to be seen by the early 
detectors: binary stars, supernova explosions, neutron 
stars, the early Universe. The estimates in this section 
are accurate only to within factors of order 1. The es- 
timates correctly show how the important observables 
scale with the properties of the systems. 

Man-made gravitational waves 

One source can be ruled out: man-made gravitational 
radiation. Imagine creating a wave generator with the 
following extreme properties. It consists of two masses 
of 10 3 kg each (a small car) at opposite ends of a beam 
10 m long. At its center the beam pivots about an 
axis. This centrifuge rotates 10 times per second. All 
the velocity is non-spherical, so v 2 N s in Equation ( |l3| ) 
is about 10 5 m 2 s -2 . The frequency of the waves will 
actually be 20 Hz, since the mass distribution of the 
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system is periodic with a period of 0.05 s, only half 
the rotation period. The wavelength of the waves will 
therefore be 1.5 x 10 7 M, about the diameter of the 
earth. In order to detect gravitational waves, not near- 
zone Newtonian gravity, the detector must be at least 
one wavelength from the source. Then the amplitude 
h can be deduced from Equation (|l3|): /i~5x 10~ 43 . 
This is far too small to contemplate detecting! 

Radiation from a spinning neutron star 

Some likely gravitational wave sources behave like the 
centrifuge, only on a grander scale. Suppose a neutron 
star of radius R spins with a frequency / and has an 
irregularity, a bump of mass m on its otherwise axially 
symmetric shape. Then the bump will emit gravita- 
tional radiation (again at frequency 2/ because it spins 
about its center of mass, so it actually has mass excesses 
on two sides of the star), and the non-spherical velocity 
will be just vjy s. — 2nRf. The radiation amplitude 
will be, from Equation (R~3|) , 



(17) 



hump ~ 2(2irRf/c) 2 Gm/rc 2 , 



and the luminosity, from Equation ( |l6| ) (assuming 
roughly 4 comparable components of Qjk contribute to 
the sum), 

L bump ~(G/5c 5 )(2irffm 2 R\ 

The radiated energy would presumably come from the 
rotational energy of the star. This would lead to a spin- 
down of the star on a timescale 



_5_ ! / Gm\ 1 rv_ 
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t spindown — TfhV / Lbump 



It is felt that neutron star crusts are not strong enough 
to support asymmetries with a mass of more than about 
m ~ 10 _5 Mq, and from this one can estimate the like- 
lihood that the observed spindown timescales of PUL- 
SARS are due to gravitational radiation. In most cases, 
it seems that gravitational wave losses cannot be the 
main spindown mechanism. 

But lower levels of radiation would still be observ- 
able by detectors under construction, and this may be 
coming from a number of stars. In particular, there 
is a class of neutron stars in X-Ray Binary Stars. 
They are accreting, and it is possible that accretion 
will create some kind of mass asymmetry or else lead to 
a rotational instability of the CFS type in the r-modes 
(see below). In either case, the stars could turn out be 
long-lived sources of gravitational waves. 

Radiation from a binary star system 

Another "centrifuge" is a binary star system. Two stars 
of the same mass M in a circular orbit of radius R 



have v 2 ^ g = GM/AR. The gravitational- wave ampli- 
tude from Equation ( |l3| ) can then be written 
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Compare this to the implications of putting Equa- 
tion @ into Equation (|). 

The gravitational-wave luminosity of such a system 
is, by a calculation analogous to that for bumps on neu- 
tron stars, 
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In this equation there appears the important constant 
c 5 /G = 3.6 x 10 52 W, a number with the dimen- 
sions of luminosity built only from fundamental con- 
stants. By comparison, the luminosity of the Sun is 
only 3.8 x 10 26 W. Close binaries can therefore radiate 
more energy in gravitational waves than in light. 

The radiation of energy by the orbital motion causes 
the orbit to shrink. The shrinking will make any ob- 
served gravitational waves increase in frequency with 
time. This is called a chirp. The timescale for this is 



''chirp 



„, 2 ,T 20GM [GM\ 

Mv /Lunary — [j^x ) 



(19) 



The binary pulsar system 
tional waves 



verifying gravita- 



This orbital shrinking has already been observed in the 
Hulse- Taylor Pulsar system, containing the radio 
pulsar PSR1913+16 and an unseen neutron star in a 
binary orbit. Discovered in 1974 by R Hulse and J 
Taylor, it has established that gravitational radiation 
is correctly described by general relativity. For their 
discovery, Hulse and Taylor received the 1993 Nobel 
Prize for Physics. 

The key to the importance of this binary system 
is that all of the important parameters of the system 
can be measured before one takes account of the orbital 
shrinking. This is because a number of post-Newtonian 
effects on the arrival time of pulses at the Earth, such 
as the precession of the position of the periastron and 
the time-dependent gravitational redshift of the pulsar 
period as it approaches and recedes from its companion, 
are measured in this system. They fully determine the 
masses and separation of the stars and the inclination 
and eccentricity of their orbit. From these numbers, 
without any free parameters, it is possible to compute 
the shrinking timescale predicted by general relativity. 
The observed rate matches the predicted rate to within 
the observational errors of less than 1%. 

The stars are in an eccentric orbit (e = 0.615) and 
both have masses of 1.4M Q . The orbital semimajor 
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axis is about 7 x 10 8 m. Equation (|T^) assumes a circu- 
lar orbit and gives a shrinking timescale of 6 x 10 10 y. 
This is an overestimate, however, partly because it is 
in any case a rough approximation, and partly because 
the timescale is very sensitive to eccentricity. With the 
observed eccentricity, a careful calculation shows that 
the expected shrinking timescale is around 4 x 10 8 y, 
consistent with observations. 



Chirping binaries 

For a circular equal-mass binary, the orbital shrinking 
timescale and the frequency of the orbit determine both 
M and R. If in addition a gravitational wave detector 
measures the wave amplitude hbi nary , then the distance 
r to the binary system can be determined. 

Remarkably, this conclusion holds even for bina- 
ries with unequal masses. In such a case, the measur- 
able mass is the chirp mass of the binary, defined as 
M := /i 3 / 5 M^ /5 , where \x is the reduced mass of the bi- 
nary system and Mt its total mass. Then the distance 
r is still measurable from the chirp rate, frequency, and 
amplitude. In other words, a chirping binary is a stan- 
dard candle in astronomy. Post-Newtonian corrections 
to the orbit, if observed in the waveform, can determine 
the individual masses of the stars and even their spins. 

Recognizing weak signals 

For ground-based detectors, all expected signals have 
amplitudes that are close to or even below the instru- 
mental noise level in the detector output. Such signals 
can nevertheless be detected with confidence if their 
waveform matches an expected waveform. The pattern 
recognition technique that will be used by detector sci- 
entists is called matched filtering. 

Matched filtering works by multiplying the output of 
the detector by a function of time (called the template) 
that represents an expected waveform, and summing 
(integrating) the result. If there is a signal matching 
the waveform buried in the noise then the output of the 
filter will be higher than expected for pure noise. 

A simple example of such a filter is the Fourier trans- 
form, which is a matched filter for a constant-frequency 
signal. The noise power in the data stream is spread out 
over the spectrum, while the power in the signal is con- 
centrated in a single frequency. This makes the signal 
easier to recognize. The improvement of the signal-to- 
noise ratio for the amplitude of the signal is propor- 
tional to the square-root of the number of cycles of the 
wave contained in the data. This is well-known for the 
Fourier transform, and it is generally true for matched 
filtering. 

Matched filtering can make big demands on com- 
putation, for several reasons. First, the arrival time of 
a short-duration signal is generally not known, so the 



template has to be multiplied into the data stream at 
each distinguishable arrival time. This is then a cor- 
relation of the template with the data stream. Nor- 
mally this is done efficiently using Fast Fourier Trans- 
form methods. 

Second, the expected signal usually depends on a 
number of unknown parameters. For example, the ra- 
diation from a binary system depends on the chirp mass 
M. , and it might arrive with an arbitrary phase. There- 
fore, many related templates must be separately applied 
to the data to cover the whole family of signals. 

Third, matched filtering enhances the signal only if 
the template stays in phase with the signal for the whole 
data set. If they go out of phase, the method begins to 
reduce the signal-to-noise ratio. For long-duration sig- 
nals, such as for low-mass neutron-star coalescing bina- 
ries or continuous-wave signals from neutron stars (see 
below), this requires the analysis of large data sets, and 
often forces the introduction of additional parameters 
to allow for small effects that can make the signal drift 
out of phase with the template. It also means that the 
method works well only if there is a good prediction of 
the form of the signal. 

Because the first signals will be weak, matched fil- 
tering will be used wherever possible. As a simple rule 
of thumb, the detectability of a signal depends on its 
effective amplitude h e ff, defined as 



h e ff = h^J N 'cycles, 



(20) 



where N cyc i es is the number of cycles in the waveform 
that are matched by the template. 

For example, the effective amplitude of the radia- 
tion from a bump on a neutron star (Equation (|l7|))will 
be ht >ump y / 2fT bs, where T a i, s is the observation time. 
In order to detect this radiation, detectors may need 
to observe for long periods, say 4 months, during which 
they accumulate billions of cycles of the waveform. Dur- 
ing this time, the star may spin down by a detectable 
amount, and the motion of the Earth introduces large 
changes in the apparent frequency of the signal, so 
matched filtering needs to be done with care and preci- 
sion. 

Another example is a binary system followed to co- 
alescence, i.e. where the chirp time in Equation ( |l9| ) is 
less than the observing time. For neutron-star binaries 
observed by ground-based detectors this will always be 
the case (see the next section), so the effective ampli- 
tude is roughly 



hchirp ^ hbin \/ f gwt chi 



GM ( GM \ 
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(21) 



where for f gw one must use twice the orbital frequency 
\J GM/R 3 /4ir. This may seem a puzzling result, be- 
cause it says that the effective amplitude of the signal 
gets smaller as the stars get closer. But this just means 
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that the signal will be more detectable if it is picked up 
earlier, since Equation (^l|) assumes that the signal is 
followed right to coalescence. If one picks up the sig- 
nal at earlier times, then there are more cycles of the 
waveform to filter for, and this naturally gives a bet- 
ter signal-to-noise ratio. This gives an advantage to 
detectors that can operate at lower frequencies. This 
has been an important consideration in the design of 
modern detectors. (See Gravitational Radiation 
Detection on Earth and in Space.) 

In general, the sensitivity of detectors will be lim- 
ited not just by detector technology, but also by the 
duration of the observation, the quality of the signal 
predictions, and even by the availability of computer 
processing power for the data analysis. 

Astronomical sources of gravitational 

waves 

Estimating the frequency 

The signals for which the best waveform predictions are 
available have narrowly defined frequencies. In some 
cases the frequency is dominated by an existing motion, 
such as the spin of a pulsar. But in most cases the 
frequency will be related to the natural frequency for a 
self-gravitating body, defined as 

fa = y/Gp/An, (22) 

where p is the mean density of mass-energy in the 
source. This is of the same order as the binary orbital 
frequency and the fundamental pulsation frequency of 
the body. 

The frequency is determined by the size R and mass 
M of the source, taking p — 3M/47ri? 3 . For a neutron 
star of mass IAMq and radius 10 km, the natural fre- 
quency is /o = 1.9 kHz. For a black hole of mass 10 Mq 
and radius 2GM/c 2 = 30 km, it is / = 1 kHz. And for 
a large black hole of mass 2.5 x 10 6 M Q , such as the one 
at the center of our Galaxy, this goes down in inverse 
proportion to the mass to fo = 4 mHz. 

Fig. ^ shows the mass-radius diagram for likely 
sources of gravitational waves. Three lines of constant 
natural frequency are plotted: fo = 10 4 Hz, /o = 1 Hz, 
and fo = 10~ 4 Hz. These are interesting frequencies 
from the point of view of observing techniques: grav- 
itational waves between 1 and 10 4 Hz are accessible 
to ground-based detectors, while lower frequencies are 
observable only from space. (See the article Gravi- 
tational Radiaton Detection on Earth and in 
Space.) Also shown is the line marking the black- hole 
boundary. This has the equation R — 2GM/c 2 . There 
are no objects below this line. This line cuts through 
the ground-based frequency band in such a way as to 
restrict ground-based instruments to looking at stellar- 
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Figure 2: Mass-radius plot for gravitational wave 
sources. 

mass objects. Nothing over a mass of about 10 4 M Q can 
radiate above 1 Hz. 

A number of typical relativistic objects arc placed 
in the diagram: a neutron star, a binary pair of neu- 
tron stars that spirals together as they orbit, some black 
holes. Two other interesting lines are drawn. The lower 
(dashed) line is the 1-year coalescence line, where the 
orbital shrinking timescale in Equation ( |l9| ) is less than 
one year. The upper (solid) line is the 1-year chirp line: 
if a binary lies below this line then its orbit will shrink 
enough to make its orbital frequency increase by a mea- 
surable amount in one year. (In a one-year observation 
one can in principle measure changes in frequency of 
1 yr-\ or 3 x 10~ 8 Hz.) 

It is clear from the figure that any binary system 
that is observed from the ground will coalesce within an 
observing time of one year. Since pulsar statistics sug- 
gest that this happens less than once every 10 5 years 
in our Galaxy, ground-based detectors must be able to 
register these events in a volume of space containing 
at least 10 6 galaxies in order to have a hope of see- 
ing occasional coalescences. When detectors reach this 
sensitivity (sometime in the first decade of the 21st cen- 
tury), then astronomers will be able to use the observed 
chirping binaries as standard candles to measure dis- 
tance scales in the Universe. 

Radiation from neutron-star normal modes 

In Fig. U there is a dot for the typical neutron-star. The 
corresponding frequency is the fundamental vibrational 
frequency of such an object. In fact, neutron stars have 
a rich spectrum of non-radial normal modes, which fall 
into several families: f-, g-, p-, w-, and r-modes have all 
been studied. If their gravitational wave emissions can 
be detected, then the details of their spectra would be 
a sensitive probe of their structure and of the Equa- 
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tion of State of Neutron Stars, in much the same 
way that Helioseismology probes the interior of the 
Sun. This is a challenge to ground-based detectors, 
which cannot yet make sensitive observations as high 
as 10 kHz. 

Radiation from gravitational collapse 

The event that forms a neutron star is the gravitational 
collapse that produces in a supernova. It is difficult 
to predict the waveform or amplitude expected from 
this event, because we have no observational evidence 
about how nonspherical the collapse event might be in 
a typical supernova: the collapse is hidden deep within 
the star. So we can only guess. For example, a gravi- 
tational wave burst might be broad-band, centered on 
1 kHz, or it might be a few cycles of radiation at a fre- 
quency anywhere between 100 Hz and 10 kHz, chirping 
up or down. The amplitude could be large, in which a 
good fraction of the energy released by the collapse is 
radiated in gravitational waves, or it could be negligi- 
bly small. It is indeed ironic that, although detecting 
supernovae was the initial goal of detector development 
when it started 4 decades ago, little more is known to- 
day about what to expect than scientists knew then. 

Radiation from r-modes 

Hot neutron stars that rotate faster than about 100- 
200 Hz appear to be unstable to the emission of gravita- 
tional radiation through amplification of their r-modes 
by the CFS mechanism. In stars colder than about 
10 s K, viscosity may be strong enough to damp out this 
instability. This instability may explain why only old, 
recycled, cold pulsars are seen at higher rotation rates. 
It also suggests that the formation of a rapidly rotating 
neutron star may be followed by a period of steady grav- 
itational radiation as the star emits angular momentum 
and spins down to its stability limit. If as few as 10% 
of all the neutron stars formed since Star Formation 
began (at a redshift of perhaps 4) went through such 
a spindown, then they may have produced a detectable 
random background of gravitational radiation. 

Interestingly, the r-modes are disturbances primar- 
ily of the fluid velocity; they have little density pertur- 
bation. Their name comes from their similarity to the 
itossby waves of oceanography. The gravitational ra- 
diation they emit is not primarily mass-quadrupole (as 
in Equation (E2)), but rather mass-current-quadrupole, 
the analog of magnetic quadrupole radiation in elec- 
tromagnetism. This is the wave counterpart of what 
is called gravitomagnetism, which is responsible for the 
Lense- Thirring effect: an extra precession of a spinning 
gyroscope as it orbits a rotating body like the Earth 
caused by the spin-spin gravitational coupling of the 
gyroscope to the Earth. 



Black holes and gravitational waves 

Black holes are regions of spacetime within which every- 
thing is trapped: light cannot escape, nor can anything 
else that moves slower than light. The boundary of 
this region is called the event horizon. This boundary 
is a dynamical surface. If any mass-energy falls into 
the hole, the area of the horizon increases. In addition, 
the horizon will generally wobble when this happens. 
These wobbles settle down quickly, emitting gravita- 
tional waves, and leaving a smooth (and slightly larger) 
horizon afterwards. 

Undisturbed black holes are time-independent and 
smooth. In fact, according to general relativity the ex- 
ternal gravitational field of such a black hole and the 
size and shape of its horizon are fully determined by 
only three numbers: the total mass, electric charge 
and angular momentum of the black hole. This black- 
hole uniqueness theorem is remarkable, considering how 
much variety there can be in the material that collapsed 
to form the black hole and that may have subsequently 
fallen in. 

Observations of the gravitational waves emitted by a 
wobbling horizon or by a particle in orbit around a black 
hole have the potential to test the uniqueness theorem 
and thereby to verify the predictions of general relativ- 
ity about the strongest possible gravitational fields. 

Astronomers now recognize that there is an abun- 
dance of black holes in the universe. Observations of 
various kinds have located black holes in X-ray binary 
systems in the Galaxy and in the centers of galaxies. 

These two classes of black holes have very differ- 
ent masses. Stellar black holes typically have masses 
of around 10M Q , and are thought to have been formed 
by the gravitational collapse of the center of a large, 
evolved Red Giant Star, perhaps in a supernova ex- 
plosion. Massive black holes in galactic centers seem to 
have masses between 10 6 and 10 10 M Q , but their history 
and method of formation are not yet understood. 

Both kinds of black hole can radiate gravitational 
waves. According to Fig. |2|, stellar black hole radia- 
tion will be in the ground-based frequency range, while 
galactic holes are detectable only from space. The ra- 
diation from a black hole typically is strongly damped, 
lasting only a few cycles about the frequency, which for 
a spherical black hole is given by Equation ( p2[ ) with 
R = 2GM/c 2 : 



Stellar-mass black holes 

Radiation from stellar black holes is expected mainly 
from coalescing binary systems, when one or both of 
the components is a black hole. Although such systems 
are thought to be rarer than systems of two neutron 
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stars, the larger mass of the black hole makes the system 
visible from a greater distance. By measuring the chirp 
mass (as discussed above) observers will recognize that 
they have a black-hole system. It is very possible that 
the first observations of binaries by interferometers will 
be of black holes. 

When a two-black-hole binary coalesces, there 
should be a burst of gravitational radiation that will 
depend in detail on the masses and spins of the objects. 
Numerical simulations of such events will be needed to 
interpret this signal, and possibly even to extract it from 
the instrumental noise of the detector. The research 
field of numerical relativity is making rapid progress, 
and it can be expected to produce informative simula- 
tions in the first few years of the 21st century, using the 
largest and fastest computers available at that time. 

Massive and supermassive black holes 

Gravitational radiation is expected from supermassive 
black holes in two ways. In one scenario, two massive 
black holes spiral together in a much more powerful 
version of the coalescence we have just discussed. The 
frequency is much lower, but the amplitude is higher. 
Equation (|2l| ) implies that the effective signal ampli- 
tude is almost linear in the masses of the holes, so that 
a signal from two 1O 6 M0 black holes will have an ampli- 
tude 10 5 times bigger than the signal from two IOMq 
holes at the same distance. Even allowing for differ- 
ences in technology, space-based detectors will be able 
to study such events with a very high signal-to-noise 
ratio no matter where in the universe they occur. 

Observations of coalescing massive black-hole bina- 
ries will therefore provide strong tests of the validity of 
general relativity in the regime of strong gravitational 
fields, provided that numerical simulations can match 
the accuracy of the observations by that time. 

The event rate for such coalescences is not easy to 
predict: it could be zero, but it may be large. It seems 
that the central core of most galaxies may contain a 
black hole of at least lO 6 Af . This is known to be true 
for our galaxy and for a number of others nearby. Su- 
permassive black holes (up to a few times 10 9 M Q ) are 
believed to power Quasi-Stellar Objects and Ac- 
tive Galaxies. There is some evidence that the mass 
of the central black hole is proportional to the mass of 
the core of the host galaxy. 

If black holes are formed with their galaxies, in a sin- 
gle spherical gravitational collapse event, and if noth- 
ing happens to them after that, then coalescences will 
never be seen. But it is believed that Galaxy Forma- 
tion probably occurred through the merger of smaller 
units, sub-galaxies of masses upwards of 10 6 M Q . If 
these units had their own black holes, then the mergers 
would have resulted in the coalescence of many of the 
holes on a timescale shorter than the present age of the 



universe. This would give an event rate of several per 
year. If the supermassive black holes were formed from 
smaller holes in a hierarchical merger scenario, then the 
event rate could be hundreds or thousands per year. It 
is likely that only space-based observations of gravita- 
tional waves will answer these questions. 

A second scenario for the production of radiation 
by massive black holes is the swallowing of a stellar- 
mass black hole or a neutron star by the large hole. 
Massive black holes exist in the middle of dense star 
clusters. The tidal disruption of main-sequence or gi- 
ant stars that stray too close to the hole is thought to 
provide the gas that powers the quasar phenomenon. 
These clusters will also contain a good number of neu- 
tron stars and stellar-mass black holes. They are too 
compact to be disrupted by the hole even if they fall 
directly into it. 

Such captures therefore emit a gravitational wave 
signal that may be approximated by studying the mo- 
tion of a "point mass" near a black hole. It will again 
emit a chirp of radiation, but in this case the orbit may 
be very eccentric. The details of the waveform encode 
information about the geometry of space-time near the 
hole. In particular, it may be possible to measure the 
mass and spin of the hole and thereby to test the unique- 
ness theorem for black holes. The event rate is not 
very dependent on the details of galaxy formation, and 
is probably high enough for many detections per year 
from a space-based detector. 

Gravitational waves from the Big Bang 

Gravitational waves have traveled almost unimpeded 
through the universe since they were generated at times 
as early as 10~ 24 s after the Big Bang. Observing them 
would provide important constraints on theories of In- 
flation and high-energy physics. 

Inflation is an attractive scenario for the early uni- 
verse because it makes the large-scale homogeneity of 
the universe easier to understand. It also provides 
a mechanism for producing initial density perturba- 
tions large enough to evolve into galaxies as the uni- 
verse expands. These perturbations are accompanied 
by gravitational-field perturbations that travel through 
the universe, redshifting in the same way that photons 
do. Today these perturbations should form a random 
background of gravitational radiation. 

The perturbations arise by parametric amplification 
of quantum fluctuations in the gravitational wave field 
that existed before inflation began. The huge expansion 
associated with inflation puts energy into these fluctu- 
ations, converting them into real gravitational waves 
with classical amplitudes. 

If inflation did not occur, then the perturbations 
that led to galaxies must have arisen in some other way, 
and it is possible that this alternative mechanism also 
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produced gravitational waves. One candidate is cos- 
mic defects, including cosmic strings and cosmic tex- 
ture. (Sec Topological Defects in Cosmology.) 
Although observations at present seem to rule cosmic 
defects out as a candidate for galaxy formation, cos- 
mic strings may nevertheless have produced observable 
gravitational waves. 

If inflation did not occur, there could also be a ther- 
mal background of gravitational waves at a tempera- 
ture similar to that of the cosmological microwave back- 
ground, but this radiation would have such a high fre- 
quency that it would not be detectable by any known 
or proposed technique. 

The random background will be detectable as a noise 
in the detector that competes with instrumental noise. 
In a single detector, such as the first space-based de- 
tector, this noise must be larger than the instrumental 
noise to be detected, and one must have great confi- 
dence in the detector in order to claim that the observed 
noise is external. This is how the cosmic microwave 
background was originally discovered in a radio tele- 
scope. 

If there are two detectors, then one can multiply 
the outputs of the two detectors together and sum (in- 
tegrate). In this way, the random wave field in one de- 
tector acts like a matched filtering template, matching 
the random field in the other detector. This allows the 
detection of noise that is below the instrumental noise 
of the individual detectors. For this to work, the two 
detectors must be close enough together to experience 
the same random wave field. In practice, the sensitivity 
of this method falls off rapidly with separation if the 
detectors are more than a wavelength apart. 

Measure of the strength of random gravitational 

waves 

When describing the strength of a random wave field, 
it is not appropriate to measure the amplitude of any 
single component. Rather, the r.m.s. amplitude of the 
field is the observable quantity. It is common to use 
an equivalent measure, the energy density p gw (f) m 
the radiation field as a function of frequency /. For a 
cosmological field, what is relevant is to normalize this 
energy density to the critical density p c required to close 
the universe. It is thus conventional to define 



This is roughly the fraction of the closure energy density 
in random gravitational waves between the frequency / 
and 2.718/. 

Current and planned detectors may reach a sensi- 
tivity of Q gw ~ 1CT 9 at 1 mHz and 1(T 10 at 40 Hz, 
but there is a possibility that backgrounds due to other 



sources (binary white dwarf systems and r-mode spin- 
down, as discussed above) could obscure a cosmological 
background at these levels. 

Predicted spectrum of cosmological radiation 

The simplest models of inflation suggest that the spec- 
trum of the gravitational wave background should be 
flat, so that fl gw is independent of frequency over a 
very large range of frequencies. In this case, the ob- 
served fluctuations in the cosmic microwave background 
radiation set a limit on gravitational radiation at ultra- 
low frequencies, and this constrains the energy density 
in the observable range (0.1 mHz to 10 kHz) to below 
about 10~ 13 of closure. This will be too small to be seen 
by the current and planned detectors on the ground or 
in space. 

But there is a great deal of room in these models for 
other spectra. The period before inflation may produce 
initial conditions for the phase of parametric amplifi- 
cation that give large amounts of radiation in the ob- 
servable frequency range. One family of models based 
in superstring theory has a spectrum that rises at high 
frequencies. If a cosmological background from inflation 
or from cosmic defects can be observed, it will contain 
important clues to the nature of the theory that unifies 
gravitation with the rest of quantum physics. 

Conclusions 

The first few years of the 21st century should see the 
first direct detections of gravitational radiation and the 
opening of the field of gravitational wave astronomy. 
Beyond that, over a period of a decade or more, one 
may expect observations to yield important and use- 
ful information about binary systems, stellar evolution, 
neutron stars, black holes, strong gravitational fields, 
and cosmology. 

If gravitational wave astronomy follows the example 
of other fields, like X-Ray Astronomy and Radio 
Astronomy, then at some level of sensitivity it will 
begin to discover sources that were completely unex- 
pected. Many scientists think the chance of this hap- 
pening early is very good, since the processes that pro- 
duce gravitational waves are so different from those that 
produce the electromagnetic radiation on which most 
present knowledge of the universe is based, and since 
more than 90% of the matter in the universe is dark 
and interacts with visible matter only through gravita- 
tion. 

Present and planned detectors are known not to be 
ideal for some kinds of gravitational wave sources. Sen- 
sitive measurements of a cosmological background of 
radiation from the big bang may not be possible with 
these instruments if the spectrum follows the predic- 
tions of "standard" inflation theory. Most of the nor- 
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mal mode oscillations of neutron stars will be very hard 
to detect, because the radiation is weak and at a high 
frequency, but the science there is compelling: neutron- 
star seismology may be the only way to probe the in- 
teriors of neutron stars and understand these complex 
and fascinating objects. Detector technology will con- 
tinually improve, and these sources provide important 
long-term goals for this field. There will clearly be much 
to do after the first observations are successfully made. 
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